clear;clc;
N = 512;
step = 1000;    

matrix = load(sprintf('./Rho_%d.dat',step));
    for j = 1:1:N
        for i = 1:1:N
            k = (j-1)*N+i;
            rho(i,j) = matrix(k,1);
        end 
    end